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Abstract 

We study the Loschmidt echo (LE) of a coupled system consisting of a central spin and its 
surrounding environment described by a general XY spin-chain model. The quantum dynamics 
of the LE is shown to be remarkably influenced by the quantum criticality of the spin chain. In 
particular, the decaying behavior of the LE is found to be controlled by the anisotropy parameter 
of the spin chain. Furthermore, we show that due to the coupling to the spin chain, the ground- 
state Berry phase for the central spin becomes nonanalytical and its derivative with respect to 
the magnetic parameter A in spin chain diverges along the critical line A = 1, which suggests an 
alternative measurement of the quantum criticality of the spin chain. 



PACS numbers: 75.10.Pq, 03.65. Vf, 05.30.Pr, 42.50.Vk 



Quantum phase transition (QPT), which is closely associated with the occurrence of non- 
analyticity of the ground-state energy as a function of the coupling parameters in the sys- 
tem's Hamiltonian[l|, are of extensive current interest, mainly in condensed matter physics 
because they are not only at the origin of unusual finite temperature properties but also 
promote the formation of new states of matter like unconventional superconductivity in 
heavy-fermion system [2|]. In the parameter space, the points of nonanalyticity of the ground- 
state energy density are referred to as critical points and define QPT. At these points one 
typically witnesses the divergence of the length associated with the two-point correlation 
function of some relevant quantum field. In experiments QPT has been extensively studied 
in the heavy-fermion compounds^, 0]. Recently, QPT has drawn a considerable interest 
in other fields of physics. More specifically QPT has been studied by analyzing pealing, 
asymptotical behavior and extremal points of various entanglement measures 3, 



The connection between geometric Berry phase (BP) and QPT for the case of spin-XY 

nnn 

model has also been studied [HJ, LLU, II 21] . through which a remarkable relation between the 
BP and criticality of spin chains is established. In addition, a characterization of QPT in 
terms of the overlap between two ground states obtained for two different values of external 



parameters has been presented 



1$. 



Another way to study quantum criticality is to investigate quantum dynamics of the 
many-body systems. Recently, Sengupta et a/.Q] have studied time evolution of the Ising 
order correlations under a time-dependent transverse field and shown that the order param- 
eter is best enhanced in the vicinity of the quantum critical point. Quan et a/.|l£| have 
studied transition dynamics of a quantum two-level system from a pure state to a mixed 
one induced by the quantum criticality of the surrounding many-body system. They have 
shown that the d ecay ing behavior of the LE is best enhanced by the QPT of the surrounding 
system. Yi et a/.[l6( have reported the relation between the Hahn spin echo of a spin-1/2 
particle and QPT in a spin chain which is coupled to the particle. It is expected that further 
work associated with the dynamical measurement of QPT via a coupling to the central probe 
system will be reported afterwards. From this aspect a thorough theoretical investigation 
of the quantum dynamics in QPT regime, including the various kinds of spin-chain models, 
is necessary and will be helpful for future experimental references. 

In this paper, we present a theoretical study of the behavior of the Loschmidt echo (LE) 
of a coupled spin system which consists of two quantum subsystems. One subsystem is char- 
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acterized by a spin- 1/2 Hamiltonian, which denotes the general two-level particles. We call 
this subsystem the central spin, in the sense that this spin plays the role of measuring appa- 
ratus. Whereas the other subsystem plays the role of surrounding many-body environment 
and is modeled by a general XY spin chain in a transverse magnetic field. The present study 
is directly motivated by the recent theoretical report Q| that the quantum critical behavior 
of environmental system strongly affects its capability of enhancing the decay of LE. Here we 
extend the Ising model used in Ref.[15] for simulating the environmental subsystem to the 
more general XY model. Compared to the Ising model, the XY model is parametrized by 
7 and A [see Eq. 1(b) below]. Two distinct critical regions appear in parameter space: the 
segment (7, A) = (0, (0, IV) for the XX spin chain and the critical line A c = 1 for the whole 
family of the XY model[l|. The behavior of decaying enhancement of the LE calculated 
in Ref . [15] can be used as a measure of the presence of the quantum criticality of the Ising 
spin chain. It remains yet to be exploited whether this decaying enhacement sustains in the 
whole critical regions for the XY model. 

The other interest in this paper is to study the BP properties of the coupled system. 
Instead of investigating the BP of the environmental XY spin chain which has been pre- 
viously studied jlo| llllll^. we focus our attention to the ground-state BP of the central 
quantum subsystem. Due to the coupling, it is expected that the quantum criticality of 
the surrounding XY spin chain will influence the BP of the central spin, which is found in 
this paper to be close proximity to the nonanalytical and divergent behavior of QPT of the 
environmental spin chain in the critical region. 

We consider a two-level quantum system (central spin) transversely coupled to a environ- 
mental spin chain which is described by the one- dimensional XY model. The corresponding 
Hamitonian is given by H = He + He + Hj, where (we take H — 1) 

H c = Ho z /2 + vo x /2, (la) 



JV 



H * = - E ( 4^*f°w + V 3 ^! + Xa! ) ' (lb) 



N 



H i = f E^f" ( lc ) 
1=1 

Here the Pauli matrices a a (a = x, y, z) and of are used to describe the central spin 



and the environmental spin-chain subsystems, respectively. The parameter A in He is the 



intensity of the magnetic filed applied along z-axis, and 7 measures the anisotropy in the 
in-plane interaction. It is well known that the XY model in Eq. (lb) encompasses two 
other well-known spin models: it turns into transverse Ising chain for 7 = 1 and the XX 
chain for 7 = 0. Hj gives the coupling between the central spin and the surrounding spin 



chain. The a 
generalization 
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2l|. 



oyed model is similar to the Hepp-Coleman model Q] or its 



As for quantum criticality in the XY model, there are two universality classes depending 
on the anisotropy 7. The critical features are characterized in terms of a critical exponent 
v defined by £ ~ | A — A c | ^ with £ representing the correlation length. For any value of 7, 
quantum criticality occurs at a critical magnetic field A c = 1. For the interval < 7 < 1 the 
model belongs to the Ising universality class characterized by the critical exponent v — 1, 
while for 7 = the model belongs to the XX universality class with v = 1/2 [lj . 

Following Ref.[15], we assume that the central spin is initially in a superposition state 
\4>s(0)) = c g \g) + c e \e), where \g) = (sin|, — cos|) T and |e) = (cos|,sin|) T with 6 = 
tan _1 (z///i) are ground and excited states of He, respectively. The coefficients c g and c e 
satisfy the normalization condition, \c g \ 2 + |c e | 2 = 1. Then the evolution of the XY spin 
chain initially prepared in \ip(0)), will split into two branches \ip a (t)) = exp(—iH a t)\ip(0)) 
(a = g, e), and the total wave function is obtained as \4>(t)) = c g \g) ® \ f g (t)) +c e \e) <S> \ip e (t)). 
Here, the evolutions of the two branch wave functions \ f a (t)) are driven, respectively, by 
the two effective Hamiltonians 

N 

H g = (g\H\g) = H E - 8^1 ~ A, (2a) 

N 

H e = (e\H\e) = H E + + A > (2b) 

1=1 

where A = ^/ fi 2 + v 2 /2 and S = gcos6/N. Obviously, both H g and H e describe the XY 
model in a transverse field, but with a tiny difference in the field strength. The central spin 
in two different states \g) and |e) will exert slightly different backactions on the surrounding 

N N 



spin chain, which manifests as two effective potentials V g = — 5^Jerf and V e = Sy^crf. This 
difference results in the decay of the LeJ^ defined as[ 



defined asfl 



L(t) = \( Va (t)\<Pe(t))\ 2 (3) 
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The LE has be en p roved to be conveniently related to depicting quantum decoherence of the 
central system^: consider the purity defined by P = r £ic{p'c) =Trc{[Tr Ep{t)] 2 } • Here 
pit) = \i/}(t))(ij}(t)\, and Tr c(e) means tracing over the degrees of freedom for the central spin 
(environmental spin chain). A straightforward calculation reveals the relationship between 
the LE and the purity as P = 1 — 2\c g c e \ 2 [1 — L(t)]|la|. This equation indicates that the 
purity depends on the initial state of the central spin and the surrounding spin chain. For 
simplicity, we assume that the spin chain subsystem begins with its ground state. In the 
following discussion, we will focus on the quantum dynamics of the LE in the different 
parameter regions. In particular, the decay problem of LE induced by the coupling of the 
central spin and its surrounding spin chain, as has been discussed in Ref.[15] for the special 
case of Ising model, will be fully studied in the (7, A)-space. 

To diagonalizc the effective Hamiltonians Hi (i = g,e), we follow the standard 



procedure^] by defining the conventional Jordan- Wigner (JW) transformation 

a l = II t 1 - 2fl m a -) i a l + 4) > ( 4a ) 



m<l 

erf = -^Hl 1 - 2a+a m ) (a; - a/~) , (4b) 

m<l 

of = 1 - 2a+a ; . (4c) 

which maps spins to one-dimensional spinless fermions with creation (annihilation) operators 
a~l (ai). After a straightforward derivation, the effective Hamiltonians read 

N 

Hi = -^2[(a^ +1 ai + a+a l+1 ) + 7(a m a, + a+a^) (5) 
l=i 

+ (A + KiS) (1 - 2a+a ; )] - «jA, 

where Kjg Kiq = 1. Next we introduce Fourier transforms of the fermionic operators 

described by d k = ^= Y,i a;e- i27r ' fc/Ar with k = — M, .., M, M = N/2. The Hamitionians (4) 
can be diagonalized by transforming the fermion operators in momentum space and then 
using the Bogoliubov transformation. The results are 

Hi = 2A fc)i (6+ 6 fc>i - 1/2) - K t A, (6) 
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where the energy spectrums (i = g,e) are given by 



A*,* = W e 2 ^ + 7 2 sin 2 ^ with e k/l = A - cos ^ + /qjj, (7) 
and the corresponding Bogoliubov-transformed fermion operators are defined by 

bk,i = cos - % sin (8) 

with angles 9 k satisfying cos9 k = e^/A^j. It is straightforward to see that the two 
sets of normal modes are related by the equation bk, e = (cos ak)bk, g — i(sin ak)b_ kg where 
= (4 e) - * } )/2. 

The ground state \G) of Hi is the vacuum of the fermionic modes described by b k ,i\G)i = 0, 

M ( d (i) 9 (i) \ 

and can be written as \G)i = Yl k =i I cos -~2~\tyk\0)-k + i sin J, where |0)fc and 

\l)k denote the vacuum and single excitation of the fcth mode, dk, respectively. Note that 
the ground state is a tensor product of states, each lying in the two-dimensional Hilbert 
space spanned by |0)/.|0)_fc and From the relationship between the two Bogoliubov 

modes bk, e and b k>g , one can see that the ground state \G) g of the effective Hamiltonian H g can 
be obtained from the ground state \G) e of H e by the transformation \G) g = rifcii( cosa fc + 
i sin a k b+ e bt kte )\G) e . 

Now we suppose that the XY spin chain is initially in the ground state of H g , i.e., 
|y(0)) = \G) g . Then our present task is to derive the explicit expression for LE. First one 
notices that the LE in Eq. (JHJ) can be rewritten as 

L(t) = |(^(t)|^(t))| 2 = |,(G|e-^|G) 9 | 2 (9) 
= \e(G\H k (cosa k -isma k b_ k ,ebk,e)e~ tHet ]J k (cosak + ismakb^ e bt k ^)\G) e \ 2 , 

where the dynamical phase in \<p g (t)) contributed by the time evolution operator e~ lHflt 
has been eliminated by the arithmetic module operation in L(t). By using the identity 
e~ lHet b k e e lHet = b k e e~ l2Ak ' et and after a straightforward derivation, one obtains the expres- 
sion for L(t) as follows 

L(t) = | e (G|nfc(cosafc-zsina fc 6_ fcie 6 fcie )(cosa fc + ie" i4Afc - e 'sina fc 6j e 6+ fcie )|G) e | 2 (10) 
= in fc (cos 2 a fc + sin 2 « fc e- i4A ^)| 2 
= Uk=i [1 - sm 2 (2a k ) sin 2 {2A k>e t)] . 
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Remarkably, the expression for L(t) based on XY spin chain is formally same as that based 
on Ising model which has been previously reported: 15]. The difference comes from the 
time-dependent phase factor, which in the present case is the energy spectrum 2Afc j6 of XY 
spin-chain characterized by the effective Hamiltonian H e , instead of Ising model given in 
Ref. [l5j . Due to the obvious difference in the energy spectrum between XY model and Ising 
model, one may expect that the behavior of the LE in the present case will include new 
features characteristic of the XY model. 

Since each factor in Eq. (|1U|) has a norm less than unity, we may expect L(t) to decrease 
to zero in the large N limit under some reasonable conditions. This kind of factorized 
structure was first discovered and systematically studied 2l| in developing the quantum 
measurement theory in classical or macroscopic limit and has been applied to analyze the 
universality of decoherence influence from environment on quantum computing)^]. Now we 
study in detail the critical behavior of of LE near the critical point A c = 1 for finite lattice 
size N of spin chain. Following Ref. , let us first make a heuristic analysis of the features 
of the LE. For a cutoff frequency K c we define the partial product for the LE 

Lc{t) = YitiFk > L c {t), (11) 

and the corresponding partial sum S{t) =lnL c = -Efc=il lnF fe |. For small k one has 

A k ^^\\-l-6\+0(k 2 ), (12) 

and 

2 , . 4vr 2 4 7 2 5 2 fc 2 , . 

sin {2ak) w w(\-s-m\ + s-ir (13) 

As a result, if K c is small enough one has 

AE(N c ) 1 2 5 2 sm 2 (2t\X-5-l\) 

where E(N C ) = An 2 N c (N c + 1)(2N C + l)/(6N 2 ). In this case, it then follows that for a fixed 
t, 

L c (t) w exp(-rt 2 ) (15) 

when A -»• A c = 1, where r = 16E(A^ C ) 7 2 5 2 /(A + 5 - l) 2 . 

From Eq. (|T5)l it may be expected that when is large enough and A is adjusted to 
the vicinity of the critical point A c = 1, the LE will exceptionally vanish with time. In the 
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FIG. 1: (color online), (a) The LE as a function of magnetic intensity A and time t for Ising 
(7 = 1.0) spin-chain size N = 100; (b) The LE as a function of time for different values of N for 
Ising spin chain. 

thermodynamic limit, i.e., the number N of sites approaching infinite while the length of 
spin chain keeping a constant, r seems to tend to zero and thus the approximate expression 
L c (t) remains unity without any decay. This implies that our heuristic analysis cannot 
apply to the case of thermodynamic limit, in which case the small-A; approximation becomes 
invalid. Thus to reveal the close relationship between the decaying behavior of LE and QPT 
which occur only in the thermodynamic limit, all fc-components of in Eq. (11) should be 
included. On the other side, for a practical system used to demonstrate the QPT-induced 
decay of the LE, the particle number N is large, but finite, and then the practical r in Eq. 
f!15|) does not vanish. 

Figure 1(a) shows the numerical result of the LE in Eq. ()10|) as a function of magnetic 
intensity A and time t for N = 100, 5 = 0.05, and 7 = 1.0 (i.e., the case of Ising model). One 
can see that when the value of A is larger or smaller than that of A c , the LE in time domain is 
characterized by an oscillatory localization behavior. When the amplitude of A approaches 
to A c , then the degree of localization of L(t) is decreased to zero. The fundamental change 
occurs at critical point of QPT, i.e., A = A c = 1. At this point, as revealed in Fig. 1(a), the 
LE evolves from unity to zero in a very short time. Figure 1(b) shows the time evolution 
of LE for different values of lattice size at critical point A = 1 of Ising model. One can see 
that the LE decays more rapidly by increasing the size N of spin chain. Also the decaying 
amplitude is increased with increasing N. 

Figure 2 shows the LE as a function of time for different values of anisotropy parameter 
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FIG. 2: The LE as a function of time for A = 1.0 and different values of anisotropy 7. The other 
parameters are chosen to be N = 100, 5 = 0.05. 

7 in the quantum critical region (A = A c ). In the extreme anisotropy limit, i.e., for the XX 
spin model (7 = 0), one can see from Fig. 2 that the LE completely remains to be unity 
during time evolution. This full localization behavior can also be seen from the analytic 
expression Eq.(14), in which r = for 7 = 0, indicating no decay in the LE, regardless of 
the variation of A and the size of spin chain. As a consequence, the purity P of the central 
spin remains unity; the coupling induced decoherence disappears for the XX spin chain. In 
this case, the quantum criticality behavior of the surrounding spin chain dose not affect the 
localization behavior of the LE for the central spin. By smoothly tuning the value of 7 little 
out of XX model, as shown in Fig. 2, the behavior of the LE begins to be characterized by an 
interplay of the decay in a short time and the oscillations in the subsequent evolution. The 
oscillations are featured by a superposition of the collapses and the revivals. The amplitude 
of the oscillations is decreased with increasing the value of 7. Further increasing the value of 
7 will, as one can see from Fig. 2, lead to the complete decay of the LE without prominent 
revivals during the whole time evolution. Therefore, the decay of the LE and its proximity 
to the quantum criticality can be tuned by the anisotropy parameter 7. 

Now we turn to study the behavior of the ground-state BP for the central spin. Due to 
the coupling, it is expected that the BP for the central spin will be profoundly influenced 
by the occurrence of QPT in spin-chain environment. 

Similar to the above discussions, it is supposed that the XY spin chain is adiabatically in 
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the ground state \G({6k})) g of H g , which is parameterized by the series {9k} in the ground 
state. Thus the effective mean-field Hamiltonian for the central spin is given by 



H eff = H s + g (G\Hj\G) 9 (16) 



(f+ff><^+r- 



In order to generate a BP for the central spin, we change the Hamiltonian by means of a 
unitary transformation: 

U(<f>) = exp ("4^) > ( 17 ) 
where is a slowly varying parameter, changing from to 2-7T. The transformed Hamiltonian 
can be written as 

H eff {4>) = C/ + (0)// e// C/(0) (18) 

f + 1 E cos ^ ) ° z + cos ^-° y sin 

The eigen-energies of the effective Hamiltonian for the central spin are given by 



E e,g - ±' 



The corresponding eigenstates are 



— cos ^e~^ J \sm l<1 



where sin# = v/2E e . 

The acquired ground-state BP for the central spin by varying from zero to 2n is given 

by 

P B = ij\\^\g) = ir(l + cose) (21) 

n \ y/\pi + Agf{\,N)]* + u*)' 

where we have defined /(A, 7, N) = Yl!k=i cos ^ ■ m the thermodynamic limit, N — > 00, 
the summation in /(A, 7, iV) can be replaced by the integral as follows: 

/(A, 7 ,N)i^ 00 =i- r A ~rr 2 • * ^ (22) 

Z7r Jo y(A — cos^J" 1 + 7^ sin </? 
10 



FIG. 3: (color online), (a) Ground-state BP of the central spin and (b) its A-derivative as a function 
of spin-chain parameters A and 7 in the thermodynamic limit. The other parameters are chosen 
to be fjL = 0.1, v = 2.0, and g = 0.5. 



The BP (3 g for the central spin is closely related with QPT of its coupled spin-chain subsys- 
tem. To manifest this, we plot in Fig. 3 the BP f3 g and its derivative d(3 g /d\ with respect 
to the field strength A as a function of spin-chain parameters A and 7. One can see that 
given the value of 7, the BP of the central spin increases with increasing the field strength 
A. After passing through the critical line A c = 1, the BP (3 g arrives at a stable value which 
turns out to be determined by the specific values of central-spin parameters fi and v. The 
nonanalytic property of BP and its A-derivative along the whole critical line can be clearly 
seen from Fig. 3. Thus a nonanalytic ground-state GP f} g and the corresponding anomaly 
in its derivative df3 g /d\ for the central spin also witness QPT of the coupling spin-chain 
subsystem. 

To help further illustration, let us consider the most discontinuous case of XX spin model 
(7 = 0). In the thermodynamic limit, the function / [Eq. (22)] occurred in the expression 
of (3 g can be obtained explicitly for 7 = as / = 1/2 — arccos(A)/7r when A < 1 and / = 1/2 
when A > 1. Thus the BP of the central spin is given by 

n ( l _|_ M+2g[l-2arccos(A)/7r] \ (j\ < 1 ) 

Jl + T J%= (A>1) 

which clearly shows a discontinuity at A = A c = 1. On the other side, one can see that the 
value of function /(A, 7, N) in j3 g is always trivial for 7 = and every finite lattice size N, 
since 9 { k g) = or 7r for every k. The difference between the finite and infinite lattice size 
can be understood, as has been first demonstrated in Ref.^|, from the two limits N — > 00 
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FIG. 4: (color online). A-dependence of ground-state BP of the central spin coupled to a XX spin 
chain (7 = 0) with different chain sizes N . The other parameters are the same as used in Fig. 3. 

and 7 — > 0. We plot in Fig. 4 the numerical results of the BP (3 g for different values of 
spin-chain size N, in comparison with the result for the thermodynamic limit. One can see 
that the BP of the central spin displays a multi-step like behavior for the small values of 
spin chain size N. By increasing N, the BP approaches towards the case of thermodynamic 
limit with nonanalyticity only at A c . We notice that the multi-step behavior of j3 g for finite 
lattice size is a unique feature of the XX model (7 = 0), and will be completely washed out 
by deviation of 7 from zero. 

To further understand the relationship between BP of the central spin and quantum 
criticality of the coupled spin chain, we calculate the derivative d/3 g /dX as a function of 
A for 7 = 1 (Ising model) and different lattice sizes. The results are plotted in Fig. 5. 
Two prominent features can be seen: (i) The derivative dj3 g /d\ of GP is peaked around 
A = 1, as in the thermodynamic limit shown in Fig. 3(b). The amplitude of the peak is 
prominently enhanced by increasing the lattice size of spin chain; (ii) The accurate position 
A m of the peak in d/3 g /d\ is changed with changing the size N of the spin chain. The 
position A m of the peak can be regarded as a pseudocritical point [24]. We show in the 
inset (red circles) in Fig. 5 the size dependence of the peak position A m for d(3 g /d\. For 
comparison, we also plot in this inset the size dependence of the peak position in A-space 
for the A-derivative of quantity /(A, 7, N). It has been shown in Ref . [1 1| that the quantity 
/(A, 7, N) is proportional to the ground-state BP for the spin chain (instead of that for the 
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FIG. 5: (color online). A-dependence of dj3 g /d\ for the central spin which is coupled to the Ising 
spin chain (7 = 1) with different chain sizes N = 15, 51, 251, 501, 00. The behavior of d/3 g /dX for 
the central spin reflects QPT of spin chain. With increasing the chain sizes, the peak becomes 
more pronounced. The inset shows the size scaling of the position of the peak occurred in d/3 g /d\ 
(circles) and in function /(A, 7, N) (squares). 

central spin discussed here) and the peak position A m in df '(A, 7 , N) / dX tends as jV~ L803 
towards the critical point. This scaling behavior of df(X, r y,N)/dX is also clearly shown in 
the inset in Fig. 5. Remarkably, compared to the scaling behavior of df(X,j,N)/dX, i.e., 
the scaling behavior of A-derivative of ground-state BP for spin chain, the peak position X m 
in df3 g /dX in the present case approaches the critical point A c more rapidly, which is verified 
by the fact that in the inset in Fig. 5 the quantity log(l — A m ) characterizing the scaling of 
d(3 g /dX curves down more rapidly than that characteristic of df{X, r y 1 N)/dX at large values 
of spin chain size N. Thus we can see that QPT of the XY spin chain is reflected faithfully 
by the behavior of the ground-state BP and its A-derivative of the coupled central spin. 

In summary, we have analyzed the behavior of the Loschmidt echo in a coupled system 
consisting of a central spin and its surrounding environment characterized by a general XY 
spin chain. The exact expression of the LE has been obtained. The relation between the 
behavior of the LE and the occurrence of QPT in spin chain has been illustrated. The 
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decay of LE, which is closely associated with the entanglement between the two coupled 
subsystems, has been shown to be monotonically modulated by the anisotropic parameter 
7 of the spin chain. At 7 = (XX model), in particular, both the heuristic analysis and 
the numerical calculation show that the LE is completely localized to be unity without any 
decay. Furthermore, we have investigated the behavior of the ground-state BP (3 g of the 
central spin. It has been shown that the behavior of (3 g and its derivative with respect to 
the magnetic intensity A of the spin chain has a direct connection with QPT of the spin- 
chain subsystem. This connection is verified by the common feature that both BP (and 
its A-derivative) of the central spin and QPT of the coupling spin chain is characterized by 
nonanalytic behavior around the critical point (or critical line) A = A c . Thus the QPT of 
the spin chain can be revealed by studying the BP behavior of the coupled central spin. 
This work was supported by CNSF No. 10544004 and 10604010. 
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